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A search for inverse beta decay electron antineutrinos has been carried out using the 825 days 
sample of solar data obtained at SK. The absence of a significant signal, that is, contributions to 
the total SK background and their angular variations has set upper bounds on a) the absolute 
flux of solar antineutrinos originated from neutrinos $-(*-B) =< 1.8 x 10^ cm~^ which is 
equivalent to an averaged conversion probability bound of P < 3.5% (SSM-BP98 model) and b) 
their differential energy spectrum, the conversion probability is smaller than 8% for all -Ee.uis > 6.5 
MeV going down the 5% level above Ee.vis ~ 10 MeV. It is shown that an antineutrino flux would 
have the net effect of enhancing the SK signal at hep neutrino energies. The magnitude of this 
enhancement would highly depend on the, otherwise rather uncertain at this moment, steepness of 
the solar neutrino spectrum at these energies. 



Introduction. The combined action of spin flavor 
precession in a magnetic field and ordinary neutrino mat- 
ter oscillations can produce an observable flux of V^ji's 
from the Sun in the case of the neutrino being a Majo- 
rana particle. In the simplest model, where a thin layer 
of highly chaotic of magnetic field is assumed at the bot- 
tom of the convective zone (i? ~ 0.7Rq), the antineu- 
trino appearance probability at the exit of the layer can 
be written as (l) (see also Refs. 



(1) 



where the parameter 1 — 2^ ~ exp(— 4ri^Ar) includes 
the layer width( Ar ^ 0.1-Rq) and the r.m.s strength of 
the field. The antineutrino flux could be large if, i.e., 
the neutrino have passed through a MSW resonance be- 
fore arriving to the layer. The MSW resonance converts 
practically all the initial 77g flux into Fu- The field fi- 



pointed out in Ref. |^ and we show in detail in this work, 
the angular distribution is of a more complicated nature 
than was naively assumed in Ref. Nevertheless, as we 
show in the detailed analysis presented here, meaning- 
ful bounds on solar antineutrino fluxes and appearance 
probability can be obtained from the SK data. 

Solar antineutrinos could have been already detected 
at SK: due to the anomalous, forward-peaked at the 
higher energies, angle distribution the excess of high en- 
ergy solar neutrino events observed in SK above 13 MeV, 
the hep spectrum anomaly could be explained, al- 

beit partially, as an excess of positrons coming from in- 
verse beta decay appearing with a direction which fakes 
that one of real solar neutrino electrons. 

Antineutrino Cross sections. In the limit of in- 
finite nucleon mass, the cross section for the reaction 
Ve+p ^ n + e+ is given by ||l^,|| a{E^) = ciEePe- The 



nally converts them into Ve. A fraction of the Ve wiU be transition matrix element ci = 2t:^ log2/m°/ti/2 is ex- 
reconverted into by mass oscillations but this recon- pressed in terms of the free neutron decay /ti/2 value, 
version is limited in this case by the chaotic character of 
the process. 

Water Cerenkov detectors such as Kamiokande and Su- 
perKamiokande (SK) which are sensitive to the Vf. — e 
elastic interaction are also capable of detecting these 
T^e's coming from the sun. Forthcoming solar neutrino 
experiments, such as SNO and Borexino are also ex- 
pected to have a high sensitivity to them. The specific 
signature of electron antineutrinos in proton containing 
materials is the inverse beta decay process: Fe -f p — > 
n -|- e^, which produces almost isotropical monoener- 
getic positrons with a relatively high cross section. An- 
tineutrino events would contribute to the background to 
forward-peaked neutrino solar events. 

The residual angular correlation between the direction 
of the incident neutrino and the resulting electron was 
proposed already in Ref. ||^ as a way of detecting so- 
lar neutrinos. More recently, it was again suggested ||] 
the same procedure to separate statistically antineutrino 
events in the SK sample. In practice however, as it was 



where the phase-space factor / = 1.71465 ± 0.00015 fol- 
lows from calculation and = 614.6±1.3 sec is 
latest published value for the neutron half- life . With 
M„ , Mp and rrie being the masses of neutron, proton and 
electron, respectively and AM = M^-Mp ~ 1.293 MeV. 

The energy Ejj of the incident neutrino is related to 
that one of the positron Ef,+ by: — E^+ + AM. From 
the values above, we obtain: ci — (9.54 ± 0.02) x 10^*'' 
cm^/MeV^. Corrections to the total cross-section due to 
weak magnetism arising from the difference in the anoma- 
lous magnetic moments of the neutron and proton and 
radiative corrections including internal bremsstrahlung 
(see Ref. |l^ and references therein) are very small and 
can be neglected for our purposes. 

Direct recoil corrections to the total cross section can 
also be neglected in principle. However, recoil corrections 
are potentially important in relating the positron and 
antineutrino energies in order to evaluate the antineu- 
trino flux. The positron spectrum is not monoenergetic 
in this case and one has to integrate over the positron 
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angular distribution to obtain the positron yield. At 
the SK experimentjbecause the antineutrino flux ^{E^) 
would typically decrease quite rapidly with energy, the 
lack of adequate corrections will systematically overesti- 
mate the positron yield. For the SK case and taking into 
account the SSM-BP98 spectrum, the eff'ect decrease 
the positron yield by 2-8% at the main visible energy 
range ~ 6 — 10 MeV. The positron yield could decrease 
up 50% at hep neutrino energies, a region where incerti- 
tudes in the total and differential spectrum are of com- 
parable size or larger. Finite energy resolution smearing 
will however diminish this correction when integrating 
over large enough energy bins: in the range 6.5 — 20 MeV 
the net positron suppression is estimated to be at the 5% 
level, increasing up 20% at hep energies. 

Indirect recoil effects through the antineutrino flux 
spectrum are also present in the angular distribution as 
we discuss below. 

The angular distribution. For low antineutrino en- 
ergies, the positron angular distribution is well-described 
by the expression 



da 



dcos9 



1 — VeOi{Eu) cos9 



(2) 



where 9 is the lab angle between the antineutrino and 
positron directions, is the positron velocity. Except 
near threshold, this velocity is nearly unity and can be 
ignored. For a linear distribution as that given by Eq.(||), 
the values for the average cosine and the a coefficient are 
related by: a — S{cos9) /ve- In the limit of inflnite nu- 
cleon mass, the coefScient a = Q^"-* = (77^ — l)/(37y^4-l) « 
0.10 is independent of E^, rj = qa/ gv is the ratio between 
axial and vector couplings of the neutron. Thus the an- 
gular distribution of the positrons is weakly backwards, 
independent of the energy above the threshold region. 

The situation greatly changes when weak magnetism 
and 0(1/M) recoil corrections are kept as is shown in 
Ref. 1^. At higher energies, terms proportional to higher 
powers of cos 6* appear. The average cosine is changed 
by a noticeable amount. It changes sign for not-so-high 
energies {Eiy ~ 12 — 13 MeV), the distribution becomes 
then forward peaked. At first order 0(1/M), we have 
(£:^>w5 MeV) 



dcost 



{l- a'^^^ cos 9 + 13 cos^ 9^ 



(3) 



where the coefficients a^^\ P can be read off from Eq.(19) 
in Ref. The average cosine can be computed from 
here yielding the following expression which is an accu- 
rate approximation valid from the threshold up E^ w 100 
MeV (/2 — Hp — Hn = 3.706 is the anomalous nuclcon 
magnetic moment, 77' = 4?7/3(l — rj)) ||]: 

(cos6l)(i) = Vea°/3 + (1 + (1 + /2)?7'a°) EJM 

~ -0.034 Ve + 2.55 X 10~^E^iMeV). (4) 



The angular differential cross section (^) is shown in 
Fig.(^) for different neutrino energies. It is important 
to notice that the angular distribution is still linear to 
a high degree. We can rewrite an effective expression 
which is a very good approximation to the first order ex- 
pression in all the energy range. It can be shown, by 
obvious arguments or, if prefered, by a rigorous Least 
Square approach, that the best linear fit to the quadratic 
expression (0) is given by 



- ^ f 1 - (E) 



where the effective constant is: 

a^y^.(£;,)-3(cos0)W 



cos 9 



(5) 



(6) 



and (cos0)(^) is given by expression (^. The angular 
distribution can be completely parametrized by a single 
quantity which is a linear function of the neutrino energy. 
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FIG. 1. The angular distribution, Eq.(^, for different an- 
tineutrino energies. Above E„ ~ 13 MeV the distribution is 
peaked forwards. 

Recoil corrections to the angle positron distribution are 
important when relating positron and antineutrino ener- 
gies. In fact they become the dominant effect at the high- 
est SK energies. The observable quantity is the positron 
angular distribution as a function of the positron energy 
itseff dNe/dcos9{Ee) - da/dcos9{Ee) x <i>^[E^{cos9)]. 
with E^{cos9) = El/{1 - E^+{1 - cos9)/M). At the 
same positron energy, the difference of the energy of the 
parent antineutrino for positrons respectively emitted at 
backward and forward energies is of order 0{1/AI) but 
increasing quadratically with the positron energy: 



A£^ - -2Ee{Ee + AM)/M, 



(7) 



at Ee~l2- 15 MeV, the difference AEjj « 0.50 MeV. 
The signal difference seen between the forward and back- 
ward directions A = is related to the average cosine. 
For our angle distribution, we have A — (3/2)(cosf?) = 
a/2. We have, schematically. 



a AT (X AA^g 
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It is estimated that, for 1 MeV positron energy bins at 
Ee 13 MeV, the effective angular positron distribution 
is forward peaked with 0.4 — 0.5 (to be compared 

with values in Fig.|l|). 

Antineutrino asymmetry and high energy neu- 
trino enhancement. In the SuperKamiokande exper- 
iment, the extraction of the solar neutrino signals from 
the final sample is basically carried out using direc- 
tional correlation of the events to the direction of the 
sun, since the recoil electrons keep the directionality of 
the incoming neutrino. The distribution of the angle be- 
tween the reconstructed direction and the direction of the 
sun presents a peak in the sun direction, cos 9 = 1. The 
flat component for cos6' <« 0.4 is considered as back- 
ground, the excess above this baseline is defined as the 
solar neutrino signal. 

In order to count number of the solar neutrino events, 
a maximum likelihood method is employed. Signal and 
background probability functions {Ps,Pbg) are defined 
which depend on angles and energies, due to the fi- 
nite, energy dependent, angular resolution of the detec- 
tor. The background probability density is obtained from 
the data. In the ideal case, the function does not de- 
pend on cos^ and P^g = 1/2. In the real case, back- 
grounds, such as 7-rays, which may have strong direc- 
tionality in the detector coordinate system are eliminated 
by fitting a multi-order, energy dependent, polynomial 
to the detector coordinate distribution and and map- 
ping it back to the direction of the sun. Clearly, this 
kind of methods does not eliminate a systematic forward- 
peaked background as antineutrinos which is virtually 
un-distinguishable from the real neutrinos. The net ef- 
fect being that the real signal is systematically underes- 
timated at low energies and, to a much higher degree, 
overestimated at higher visible energies. As it is im- 
plied by the results to be presented below where we use 
the standard model neutrino spectrum, the net effect is 
too small in practice to explain the experimental excess 
10,^ . It is intriguing anyway that the vanishing and sign 
change of the antineutrino asymmetry occurs precisely 
at right position of the high energy neutrino spectrum 
« 12.5 MeV. Irrespective of the total flux, the situ- 
ation could change if the spectrum profile, through the 
quantity | d^/dEi, |, highly deviates from the assumed 
standard value. 

Experimental angular distributions. The direc- 
tion of the positron from antineutrino capture is smeared 
away by the angular resolution of the detector. This an- 
gular resolution is energy dependent and due to it, the 
effective, experimentally detectable, asymmetry parame- 
ter Qfe// is slightly smaller than the theoretical one and 
includes an additional small energy dependence. 

We have performed a simple Monte Carlo simulation 
in order to compute the influence of such effect in the 
SK data. We have found convenient to parameterize the 
SK angular resolution in terms of the Beta distribution 



e(6')m,n ~ 9"^~^{1 — 9)^~^, the parameters m, n has been 
obtained from a fit to the data presented in Ref. ||l^,|l5 . 
In particular the variance corresponding to the Beta dis- 
tribution can directly be obtained from a fit to the energy 
dependent angular resolution presented in Ref. jlSj . Full 
concrete expressions for as will be presented in Ref. [ p^ . 
We have found that the energy dependence of aeff can 
be parametrized as a simple linear expression. 

The main conclusion is that the net effect of the finite 
angular resolution is the smearing of the initial angular 
distribution consistent with a reduction of no more 15% 
on aeff at the lowest energies observable at SK. At the 
highest energies, the effect of the detection angular reso- 
lution decreases becoming negligible. 

Expected antineutrino flux. The average number 
of positrons Ni which are detected per visible energy bin 
AEi is given by the convolution of different quantities 

N.^Qo f dEe f dE:e{Ee)f{Ee,El) 

J AEi Jo 

poo 

X / dETjF{ETj)a{E-,E^J (9) 

where Qo is a normalization constant accounting for the 
fiducial volume and live time, F[Ey) is the flux of solar 
antineutrinos per unit energy at the detector. When ig- 
noring recoil effects the antineutrino capture cross section 
a{E-u,El) is simply given by a{E-,El) = a{E-)5{EY7 - 
El) with a{E-) given as before. The functions e(i?e) 
and f{Ee,El) are, respectively, the detection efficiency 
and the energy resolution function which includes trig- 
ger and VD efficiencies. Expressions for the functions 
e, / has been obtained by us using the data presented 
in Refs. |l^Jl^. Further analysis and explicit expres- 
sions can be found in Ref. The antineutrino flux, 
production probability P-p^{E^) and the solar-originated 

neutrino flux F are related by F = F x Pp^{E^). 

Results. We have analyzed the data corresponding to 
the full energy range 6.5 — 20 MeV from the first 504 days 
of operation of SK as reported in Ref. jl^ as well as the 
preliminary LE-trigger data corresponding to the first 
703 and 825 days of operation |]l6| , |l7| . In addition we 
have analyzed data corresponding to individual energy 
bins of 0.5 MeV interval as presented in Ref. [|3|. Note 
that, in the high energy end, individual energy bins 14 — 
15, 15 — 16 and 16 — 20 MeV are also considered. 

Limits from integrated background data. The results 
obtained from the analysis of the angle-integrated back- 
ground data over the full energy range are summarized 
in Table (||) (see second column). From the condition 
Ni < Backexp, a model independent upper limit on the 
flux of antineutrinos originated from neutrinos is ob- 
tained $-(®B) < 7.5 X 10''' CTO~2 s^^. Note that this limit 
would include antineutrinos with energies covering all the 
energy range. This number is equivalent to an upper 
bound P < 14% CI95% on the energy averaged neutrino- 
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antineutrino conversion probability (SSM-BP98). Note 
that much more strict bounds will be obtained below. 

The corresponding results obtained from the analysis 
of individual energy bins are shown in the Figs.(H). In 
Fig. (^) (top) we show per each positron visible energy bin 
(width ±0.25MeV, except higher energies) the observed 
flat background, the maximal SSM (BP98, see Ref. |l8)) 
expected positron signal (supposing P — 1). In the bot- 
tom figure (solid circles) we present upper limits on P 
which are inferred. We observe that the average con- 
version probability is below 8 — 10% for energies above 

9 MeV and below the 5% line for visible energies above 

10 MeV. In Fig.(^ we present the same information in a 
slightly different, smoother way, for integrated bins above 
a certain threshold: the solid circles represent upper lim- 
its on P obtained from the consideration of the global 
background. 

The limits from global background are complemented, 
specially in the lower energies, with those obtained from 
consideration of the their angle distribution as follows. 



Live Time 






504 days 
703 days 
825 days 


7.05 ± 0.03 
6.98 ±0.03 
6.87 ±0.02 


< 7.7 X lO"* 

< 7.6 X 10^ 

< 7.5 X 10^ 



TABLE I. SK live time data. BackExp{evt/kt /day): 
observed flat background. Limits on the antineutrino flux 
^—{^B) (cm~^ s~^). We have supposed for the antineu- 
trino a solar neutrino spectrum with a ®B absolute flux 



$.(»B) = 5.15l°:?« X 10« cm-^ s'^ 
would yield a maximal quantity of ■ 
trinos observable in SK. 



|18|. This absolute flux 
8.2 evt/kt/day antineu- 
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FIG. 3. Upper limits on antineutrino conversion probabil- 
ity as a function of -E^+.t?!- Fig-(i) 

Limits from fit to the angular distribution. Bounds on 
the antineutrino flux have been extracted from the fit 
to the counting yield in the angular region where "back- 
ground" events from ue interactions can be neglected. 
For this purpose the use of the linear expression given by 
Eq.(||) instead of the, only marginally more exact Eq.(||), 
is especially advantageous because allows for a trivial av- 
eraging over a finite energy range. The averaging of a 
complex energy-dependent expression is traded off by the 
averaging of a single parameter which depends linearly on 
the energy. An small solar neutrino contribution is still 
expected at cos 6* ~ 0.5. We have explicitly checked that 
the results to be presented below, which use the available 
information about the slope sign, are independent of the 
concrete angle range used for fitting, in practice we have 
found rather stable results when we fit any number of 
bins up cosine values in the range cos9 < 0.25 — 0.5. 
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FIG. 2. Top: observed flat background and maximal ex- 
pected positron signal. Bottom: upper limits on antineutrino 
conversion probability obtained from global counting (solid 
circles) and angle fit (empty circles). 



Live Time 


(■P("»oxp 


P<Po 




504 days 




1.6 X lO"-" 




0.075 




< 3.4 X lO'' 


700 days 




8.8 X 10"* 




0.040 




< 2.1 X 10^ 


825 days 




7.7 X 10~* 




0.035 




< 1.8 X 10^ 


TABLE 


II. 


(^'(a»exp: 


slope 


fit to 


the 


experimental 



angle distributions {Eth > 6.5 MeV), CI95% (unified 
Feldman-Cousins Approach with the boundary condition 
{P{a))Exp > ). Limits on the antineutrino probability sup- 
posing {oeff) = 0.022 (see text). Bounds on the antineutrino 
flux ^-CB) (cm-2 s"^) (see Table(@)). 

In order to draw a limit corresponding to a visible en- 
ergy bin of width AEi, the fitted slope b must be com- 
pared with the expected value aoxp = Qo/M {o'^j j) AEi, 

where {a^^fj)AEi is the energy average of Eq.(||) over the 
energy range AiS^ and M is the number of bins in which 
the full angle range is partitioned. From Eqs.(^) and (|^), 
we obtain 

< aj/ >A£;.-= -0.092 + 2.55 x IQ-^{E^+)ae,. 

For example, for the fuU range 6.5-20 MeV, {E^+) ~ 10.0 
MeV and the value for the observable distribution param- 
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eter {oi'^ff) — —0.022. Note that this value is far from 
the naive value for q;° = —0.1 quoted above. Note in ad- 
dition that detector angular resolution effects have been 
ignored here, according to estimations (see MC simula- 
tions above) they would induce a further, modest, smear- 
ing out of about 10 — 15% which conservatively has been 
included in the results. At higher energies and for small 
bin widths the recoil corrections are important and have 
throughly been included. 

The results obtained from the fit to the full energy 
range distributions are shown in Table (||). An impor- 
tant improvement with respect to the bounds derived 
from the global background is obtained. The model in- 
dependent upper limit on the antineutrino flux is now 
<i>p-(*i?) < 1.8 X 10^ cm~^ s~^, which corresponds to an 
upper bound P < 3.5% CI95. In order to extract CI in- 
tervals, we have used the known information about the 
sign of the expected slope: the unified Fcldman- Cousins 
approach for dealing with boundary problems |l9| ] has 
been used for this purpose. 

A similar analysis was performed for each of the indi- 
vidual energy bins and for the cumulated distributions. 
The results are shown in Figs. (|^) (bottom figure, empty 
circles) and (||) (empty circles). We see from the figures 
the complementary nature of both type of analysis: sen- 
sible limits are obtained only for energies down 10 MeV, 
the region where the bounds obtained from consideration 
of the global background are weaker. One draws a com- 
bined limit of P < 8% CI95% for the full energy range. 
The bounds are in fact lower than the 5% level for the 
great part of the interval. 

In summary, there are no indications of electron an- 
tineutrinos in the up-to-date SK data. This negative re- 
sult sets upper bounds on the total and differential so- 
lar antineutrino spectrum. We obtain an upper bound 
P^^v <~ 3.5% CI95% above the 6.5 MeV threshold. As 
a function of the energy, the stricter limits (~ 2%) are 
obtained for the highest positron visible energies. 

This work has been supported by DGICYT under 
Grant PB95-1077 and by a DGICYT-MEC contract at 
Univ. de Valencia. 
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